Suppression or enhancement of the Fulde-Ferrell-Larkin-Ovchinnikov order in a 
one-dimensional optical lattice with particle correlated tunnelling 
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We study through controlled numerical simulation the ground state properties of spin-polarized 
strongly interacting fermi gas in an anisotropic optical lattice, which is described by an effective 
one-dimensional general Hubbard model with particle correlated hopping rate. We show that the 
Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) type of state, while enhanced by a negative correlated 
hopping rate, can be completely suppressed by positive particle correlated hopping, yielding to an 
unusual magnetic phase even for particles with on-site attractive interaction We also find several 
different phase separation patterns for these atoms in an inhomogeneous harmonic trap, depending 
on the correlated hopping rate. 



PACS numbers: 

The recent experimental observation of normal- 
superfluidity transition in spin-polarized ultra-cold Fermi 
gas has triggered tremendous interest in this system 
0) S H, BL Q • In a spin-polarized gas, the competition 
between the Cooper pairing and the mismatch of Fermi 
surfaces of different spin species could lead to some ex- 
otic quantum phases @, UlU- On this regard, an exam- 
ple with particular interest is the Fulde-Ferrell-Larkin- 
Ovchinnikov (FFLO) state, which spontaneously breaks 
the translational and rotational symmetry in space @, 0] • 
The stability region of the FFLO state in three dimen- 
sions is unfortunately very narrow 0, [13] , and so far it has 
not been observed in the cold atomic gas. ft is much eas- 
ier to observe the FFLO-like state in a one-dimensional 
(ID) or quasi-lD gas due to the Fermi surface nesting 
[O, E3, US]- The FFLO type of state has been found 
in ID systems with a number of theoretical studies on 
the Hubbard model and its continuous version, including 
the bosonization approach [l4| , the Bethe ansatz solution 
combined with the mean field method [ll], [ijj , and the 
numerical simulations based on the density matrix renor- 
malization group (DMRG) 0, El E3 or the quantum 
Monte-Carlo 

To test the theoretical predictions and observe the 
FFLO-like states, one needs to have strongly interact- 
ing atomic gas confined in an anisotropic optical lattice, 
where the hopping is basically along only one spatial di- 
rection. The strongly interaction introduced by the Fes- 
hbach resonance is important for experimentally achiev- 
ing the low-temperature phase of the system. However, 
this strongly interaction causes additional complications 
to the theoretical model Hamiltonian. Because of the 
multi-band population and the atomic collision over the 
neighboring sites, the system cannot be described a con- 
ventional single-band Hubbard model any more, fnstead, 
based on the analysis of the local Hilbert space struc- 
ture and the symmetry argument, it is shown in (l9l. l2(i| 
that the effective Hamiltonian for strongly interacting 
fermions in an optical lattice is described by a general 
Hubbard model with particle correlated hopping rates. 



In this paper, we investigate the possibility to observe 
the FFLO-like state in the ID general Hubbard model 
(GHM) with population imbalance. Note that in an ex- 
actly ID configuration, the state can have only a quasi- 
long-range (QLO) order with divergence in the corre- 
sponding susceptibility. However, this QLO order can be 
stabilized to a true-long-range order in a quasi- ID config- 
uration with small transverse tunneling [2l|. So we can 
still use the leading QLO order to character the phase of 
the system. We show that the particle correlated hopping 
in the general Hubbard model has significant influence 
on the FFLO phase. While a negative particle correlated 
hopping rate enhances the stability of the FFLO state, 
a sufficiently large positive particle correlated hopping 
rate can completely suppress the FFLO phase and sta- 
bilize an unusual spin-density wave (SDW) state (with 
quasi-long range antiferromagnetic order) even for the 
atoms with on-site attractive interaction. The solution 
here is based on controlled numerical simulation with the 
tensor network algorithm in the thermodynamical limit 
[12] (a variation of the DMRG method 23]). We also 
give a analysis of the phase separation pattern for this 
system in a weak global harmonic trap. While the mean 
field approach finds a large region of the FFLO state sur- 
rounded by a fully paired superfluid state in the wing pl| , 
the exact numerical simulation shows that the parameter 
window for such a region is very narrow. For the case of 
the general Hubbard model, more intriguing phase sepa- 
ration patterns could occur. For instance, with positive 
particle correlation hopping rate, one can have a SDW 
state at the core and a fully paired superfluid state in 
the wing, while in between the state of the system shows 
some remanence of the FFLO-type order. 

For strongly interacting fermions near a wide Feshbach 
resonance, the effective Hamiltonian is given by the fol- 
lowing general Hubbard model H(| 
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where m a = a ia .ai a and = n.^ + are particle num- 
ber operators, \i a stands for the chemical potential for 
the spin-a species, denotes the neighboring sites, 

and a ia is the creation operator to generate a fermion on 
the site i with the spin index a. The symbol a stands for 
(|, |) given a = (f , J,). The Sg and St terms in the Hamil- 
tonian describe the particle correlated hopping, which 
come from the contributions of the multi-band processes 
and the direct neighboring coupling of the atoms. As one 
moves far away from the Feshbach resonance, the corre- 
lated hopping terms vanish and the Hamiltonian in Eq. 
(1) reduces to the conventional Hubbard model. How- 
ever, near the resonance, Sg and St can be significant 
compared with the atomic tunneling rate t. It is inter- 
esting to note that a similar form of the Hamiltonian 
has also been proposed for the high-Tc cuprate materi- 
als [13], so the interest in this model is not limited to 
the system of the strongly interacting atomic gas. In Eq. 
(1), we keep /Zj and different, which accounts for the 
population imbalance in the two spin components. 

In this work, we consider an anisotropic optical lat- 
tice where the hopping along the x,y directions are 
suppressed by the potential barriers. For this effective 
one-dimensional system, the ground state can be found 
through the time-evolving block decimation (TEBD) 
algorithm [2!q |. For this algorithm, first we trans- 
fer the fermions to effective spins through the Jodan- 
Wigner transformation [25]. Each site has four possi- 
ble states, denoted by \i s ) for the sth site. The coef- 
ficient Cj 1 ...j n of the ground-state wave function |^) = 



<L,i 1= l Z^i„=i L »i. 
trix product form: 
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is expressed in the ma- 
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where T^ ls denotes the matrix associated with site-s 
with the matrix dimension \- The lattice is bipartite, 
and for calculation in the thermodynamical limit, we as- 
sume a translational symmetry of the matrix for each 
sublattice [13, Hfl . The matrix is optimized through min- 
imization of the energy with the imaginary time evolu- 
tion by the Hamiltonian H . To identify different orders, 
we calculate the real space spin (S r ), density (D r ), and 
pair (P r ) correlations and their Fourier transformations 
Xk = l/v r /V/^ r=0 1 X r cos(fcr), where M is the number 
of sites involved in the transformation and X stands for 
S, D, or P correlation 25]. The real space correlation 
functions are defined by 

cm _ I m m \ / c "i\ / a "i \ 
J r — \ a i a i+r/ \ a i /\ B i+r/i 

D r = {rhim+r) - (ni)(n i+r ), (3) 

P r = (a iT a u aj +ri at +rT ) - (a iT a ii )(aT +ri at +rT ), 

where the spin operators associated with site i is given 
by s™ = a] Q 0-^0^/2 with a and (3 =1, | and a m (m = 
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FIG. 1: Fourier transform of the correlation functions for the 
general Hubbard model at half-filling. For plots (a)-(d), U = 
— 8t, while U — —2t for (e)-(h). For all the plots, the spin 
polarization is chosen to be the same at p — 0.4. The solid, 
dashed, and dash-dotted curves correspond to Sg — 0, — 0.5t, 
and t, respectively. The parameter St is taken as St — —25g. 



x,y,z) standing for the Pauli matrices. In the calculation, 
we typically take \ = 80 for the matrix dimension and 
M =100 for the Fourier transform (the results have been 
well converged with the above choice of the parameters). 

In Fig. 1, we investigate how the particle correlated 
hopping rate in the GHM influences the stability of the 
FFLO state and the phase of the system. In our cal- 
culation, we take attractive interaction with U = — 8t 
and — 2t respectively for the figures on the left and the 
right panels, and choose a fixed population imbalance 
with p = (JVf — Ni) I (Nf +Ni) = 0.4. First, for the con- 
ventional Hubbard model with Sg = St = 0, the pair 
correlation Pk peaks at nonzero momenta k with popu- 
lation imbalance, which is a signature of the FFLO like 
state. The result thus confirms the previous theoretical 
predictions. Then, we tune the particle correlated hop- 
ping rate Sg. With a negative Sg, as shown in Fig. I , 
the peaks of Pk at non-zero momenta become sharper, 
which means that the FFLO order gets enhanced (the 
FFLO is the leading quasi-long-range order in this ID 
configuration). However, with a positive Sg, the FFLO 
peaks in the pair correlation Pk get suppressed, and at 
Sg = t, the peaks in the spin correlation Sk get more 
prominent. The peaks in Sk are particularly sharp with 
a small on-site attraction U = —2t, and the spin density 
wave emerges as the leading quasi-long range order for 
the system. Although the spin density wave order is ex- 
pected for the repulsive Hubbard model, its appearance 
in the case of attractive interaction is unusual since at- 
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FIG. 2: Pair correlation functions for the conventional Hub- 
bard model at U = —8t. The difference in chemical po- 
tential for spin-up and spin-down species are chosen to be 
Sfj, = (hi — /tf)/2 = 2 At, 2.5t, and 3i for plots (a), (b), and 
(c), respectively. The solid, dashed, and dash-dotted lines in 
the plots correspond to /i = (/if + /if)/2 = — 4f (half-filling), 
— 4.3i, and —4At, respectively. The insets show the spin po- 
larization p corresponding to each curve. 



traction normally favors some pairing order. The result 
of this simulation shows that the particle correlated hop- 
ping can cause some qualitative change to the phase of 
the system. 

In the next step, we analyze the possible phase sep- 
aration patterns for the polarized fermi gas in an inho- 
mogeneous global trap V (r). When the trap V (r) is 
slowly varying from site to site, we can take the local 
density approximation where V (r) basically decreases 
the local chemical potential from \i (/i = (/if + /tf)/2) 
to fj, — V (r) as one moves from the trap center to the 
edge. So to investigate the qualitative phase separation 
pattern, we can fix the population imbalance p (and thus 
also Sfi — (/if — /if )/2) at the center, and calculate what 
kind of phases can emerge as one moves to the trap edge. 
For the conventional Hubbard model, the mean-field the- 
ory has predicted a large parameter region where one 
has a FFLO-like phase at the center, surrounded by a 
non-polarizing BCS phase at the edge [HI- This separa- 
tion pattern is somewhat unusual as in the 3-dimensional 
case, one can only has a BCS phase at the trap center 
surrounded by other phases (such as FFLO or polarized 
normal state) at the edge [^, [l(| • The recent DMRG cal- 
culation however does not find any evidence of this type 
of phase separation pattern predicted by the mean-field 
theory [15| . To resolve this problem, we have performed 
more extensive calculation over a large parameter win- 
dow. We indeed find that this type of phase separation 
is possible, however, it exists only for a very narrow pa- 
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FIG. 3: The correlation functions for the general Hubbard 
model at U/t = —2, 5g/t — 1 and St/t — —2. Difference 
in chemical potentials for spin-up and spin-down species are 
chosen to be Sfi = (/tf — fi\)l2 = 1.2t and 2At for the left 
and right column, respectively. The solid, dashed, and dash- 
dotted lines correspond to /i = (/if +/tf )/2 = —t (half- filling) , 
—2.5t, and — 3.3t, respectively. 

rameter window with the region much smaller than the 
one predicted by the mean-field theory. 

In Fig. 2, we show different kinds of phase separa- 
tion patterns for the conventional Hubbard model with 
U = —8t. We take the density about half filling at the 
trap center and vary the population imbalance there. 
With a tiny population imbalance at the center (for in- 
stance with p — 0.054 in Fig. 2a, with the corresponding 
Sfi = 2 At), when one moves to the edge, the population 
imbalance p decreases with decrease of /i (see the insert of 
Fig. 2a), and a non-polarizing BCS state emerges at the 
edge. However, when one slightly increase the popula- 
tion imbalance at the center (for instance with p = 0.093 
in Fig. 2b, with the corresponding <5/i = 2.5i), the pop- 
ulation imbalance increases with a decreasing \i as one 
moves to the edge, and the state finally goes to a po- 
larized normal phase (with no peaks in the correlation 
function, not shown in Fig. 2). This is different from the 
mean- field prediction where one could get a non- 
polarizing BCS state at the edge as long as the central 
polarization p is below 0.2. In Fig. 2c, we show the result 
when the central polarization p ss 0.2,. Clearly, the po- 
larization is increasing as one moves to the edge (similar 
to Fig. 2b), and there is no possibility of a BCS state 
there. 

In Fig. 3, we investigate the phase separation pat- 
tern for the GHM with particle correlated hopping. The 
phase separation pattern is of particular interest for the 
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case of a positive particle correlated hopping rate Sg. In 
this case, at the trap center we have a spin density wave 
state, while at the edge the state depends on the overall 
population imbalance of the system. In the case of a large 
population imbalance (the right panel of Fig. 3), we get 
a polarized normal state at the edge. However, in the 
case of a small population imbalance, the spin-density 
wave state can be accompanied by a non-polarizing BCS 
state at the edge, as evidenced by the sharp peak in the 
pair correlation P k in Fig. 3d. With a negative particle 
correlated hopping rate Sg, the phase separation pattern 
is qualitatively similar to the case of the conventional 
Hubbard model, and thus not shown in Fig. 3. 

In summary, through well controlled numerical simu- 



lation we have investigated the properties of the attrac- 
tive general Hubbard model under population imbalance, 
which describes the strongly interacting fermi gas in a 
one-dimensional optical lattice. The FFLO type of or- 
der gets either enhanced or suppressed depending on the 
sign of the particle correlated hopping rate Sg. When 
the FFLO state is fully suppressed with a positive Sg, we 
get an unusual spin density wave state. We also inves- 
tigate the phase separation pattern of the system under 
a weak inhomogeneous trap, and find several different 
phase separation patterns depending on the polarization 
of the system and the particle correlated hopping rates. 
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